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ABSTRACT 

Correlations are investigated of the outflow strength of quasars, as measured 
by the blueshift and asymmetry index (BAI) of the CiV line (Wang et al. 2011), 
with intensities and ratios of broad emission lines, based on composite quasar 
spectra built from the Sloan Digital Sky Survey. We find that most of the line 
ratios of other ions to C IV prominently increases with BAI. These behaviors can 
be well understood in the context of increasing metallicity with BAI. The strength 
of dominant coolant, CiV line, decreases and weak collisionally excited lines 
increase with gas metallicity as a result of the competition between different line 
coolants. Using SiIV+OlV]/ClV as an indicator of gas metallicity, we present, 
for the first time, a strong correlation between the metallicitiy and the outflow 
strength of quasars over a wide range of 1.7 to 6.9 times solar abundance. Our 
result implies that the metallicity plays an important role in the formation of 
quasar outflows, likely via affecting outflow acceleration. This effect may have 
a profound impact on galaxy evolution via momentum feedback and chemical 
enrichment. 

Subject headings: galaxies: abundances — galaxies: nuclei — line: formation — 
line: profiles — quasars: emission lines — quasars: general 

1. Introduction 

Outflows appear to be a common phenomenon in quasars. On small scale, mass loss 
due to outflows is an important component in the overall structure of quasars (Crenshaw et 
al. 2003, and references therein). Recently, outflows have been regarded as one of the main 
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components in the unification of various subclasses of quasars (Richards et al. 2011). On 
large scales, outflows are considered to be able to clear out inter-stellar medium (ISM; Silk & 
Rees 1998; Fabian 1999; Di Matteo et al. 2005) and enrich the inter-galactic medium (IGM; 
Friaca & Terlevich 1998), and consequently affect the star formation in the host galaxies. 
Quasar outflows are important for our understanding of the co-evolution of galaxies and 
their central black holes. 

There exist, in quasar spectra, prominent features associated with outflows, e.g. broad 
absorption lines (BALs) and blueshifted broad emission lines (BELs). Studies on the two 
different features reveal quite similar dependencies of outflow properties on some fundamental 
quantities. BAL quasars with larger outflow velocity tend to have higher Eddington ratio 
and weaker intrinsic X-ray relative to UV (Ganguly et al. 2007; Fan et al. 2009). And the 
fraction of BAL quasars also increases with Eddington ratio (Ganguly et al. 2007; Zhang 
et al. 2010). The blueshift of CiV BELs is more strongly correlated with the Eddington 
ratio than the BAL properties (Wang et al. 2011, hereafter Wangll) and also depend on the 
shape of the ionizing continuum (Leighly & Moore 2004; Richards et al. 2011). Furthermore, 
the blueshift is apparently stronger in radio-quiet quasars than in radio-loud ones (Marziani 
et al. 1996). These observational results provide strong support for the scenario that the 
outflow is directly related to the accretion process and driven by radiation pressure. 

Is the launch of outflows affected by other factors, in particular the environment where 
quasars reside in? In this letter, we uncover a strong correlation between the outflow strength 
and the quasar metallicity. This correlation is actually not unexpected considering that 
high metallicity can increase the acceleration of outflows. One good example with strong 
outflows and enriched metallicity is the ultraluminous infrared quasar Q132H-058 (Wang et 
al. 2009). Since metals are produced by star formation process, this correlation strongly 
suggests that star formation in galaxies can affect the evolution of the central black holes. 
The connection between galaxies and black holes might be more intimate and complex than 
previous expectation. 

This letter is organized as follows. In Section [2] we show the sample selection and present 
composite spectra as a function of the blueshift and asymmetry index (BAI) of the CiV line, 
which has been shown to be a good indicator of the strength of quasar outflows (Wangll). 
We also measure the line flux from the composites. In Section[3]we present strong correlations 
of various line flux ratios with BAI. These results provide strong evidence to support that 
outflow strength increases with the metallicity. Finally, we summarize the results and discuss 
their implications in Section HJ 
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2. Sample Selection and Composite Spectra 

2.1. Sample Selection 

We select quasars in the redshift range 1.7 < z < 4.0 from the Fifth Data Release of 
the Sloan Digital Sky Survey (SDSS; Schneider et al. 2007). The redshift range is chosen 
in such a way that the C IV line falls in the wavelength coverage of the SDSS spectrograph. 
To ensure reliable measurements of emission line parameters, we select objects with median 
signal-to-noise ratio (S/N)>7 per pixel in the CiV (1450-1700A) spectral region. We discard 
BAL quasars as cataloged by Scaringi et al. (2009), since BALs would significantly modify 
the Civ line profile. 12844 quasars meet these criteria. The SDSS spectra are transformed 
to the rest frame using the improved redshifts for SDSS quasars as computed by Hewett & 
Wild (2010, hereafter HWIO). 

To measure the C IV line, we first fit the local continuum with a power-law in two wave- 
length windows near 1450 A and 1695 A. After subtracting the continuum, we fit the residual 
spectrum around CiV with two Gaussians. Since the red wing of CiV is contaminated by 
Hell, only the spectral region of 1450-1580A is considered in fitting the CiV line. The fitting 
results for most of the objects are reasonable according to our visual inspection. However, a 
small fraction of objects cannot be well fitted by our automated procedure. To minimize the 
number of outliers, we eliminate objects with unacceptable fitting, i.e. x^/d-o.f> 1.5. This 
restriction leads to a final working sample of 11268 quasars. 

Following Wangll, we use the blueshift and asymmetry index to measure the deviation 
of the C IV line from an unshifted and symmetric profile. BAI is defined as the flux ratio of 
the blue part to the total profile, where the blue part is the part of C IV line at wavelengths 
short of 1549. 06A, the rest-frame wavelength of the CiV doublets. The BAI estimation may 
be affected by the accuracy of the redshifts which we use to transform the observed spectra 
to the rest frame. To reduce the possible uncertainty of BAI thus introduced, we adopt the 
redshifts provided by HWIO, which were derived by cross correlating observed spectra with 
a carefully-constructed template. It has been demonstrated by Wangll that the redshifts of 
HWIO do not introduce any significant bias in measuring BAI. 



2.2. Composite-Spectra Construction and Emission-Line Measurement 

With increased S/N, composite spectra are useful for studying weak lines and the overall 
properties of quasars by averaging out object-to-object variations (Vanden Berk et al. 2001, 
hereafter VBOl; Nagao et al. 2006, hereafter N06). Since we are interested in the relative 



-4- 



intensities of BELs, the composites are constructed using the arithmetic mean method. We 
first divide the quasar sample into five equally-sized subsamples according to their BAI. 
For each quasar, we use the HWIO redshift to deredshift the spectrum. The spectrum is 
then normalized at 1450A, rebinned into the same wavelength bins. In Figure [H we show 
the composite spectra covering the wavelength range from 1000 to 3000A. In our sample, 
BAI is weakly correlated with redshift and luminosity, so the number of quasars in each 
wavelength bin as a function of wavelength is similar for different composites (Figure [2]). 
At least 220 and 940 spectra are combined in the bins of the composites near lOOOA and 
3000A, respectively. This enables detections of weak features in the composite spectra at 
high statistical significance. Following VBOl and Laor et al. (1997), we locate and identify 
most of significant line features in the composites. 

As Civ BAI increases, other high ionization lines, such as Hell, Nv and OVI, also 
significantly shift towards shorter wavelength. Especially, the blue wing of these lines are 
prominently enhanced as BAI increases, suggesting that BAI is a good indicator of outflow 
strength. While, in our composites, low ionization lines, such as Ol, Cm] and Mgll, don't 
exhibit apparent shift with respect to the rest frame. It is consistent with previous findings 
(Marziani et al. 1996; Richards et al. 2002; Shen et al. 2007). These results support that 
there is no significant bias in HWIO redshifts and BAI. 

We then measure the fluxes of the metal lines (or line complex) listed in Tab. [1] from 
the composites. Generally, we first fit the continuum with a power-law in two continuum 
windows. After subtracting the continuum, we fit the residual line spectrum with multi- 
Gaussians model. The central wavelengths of the continuum windows and the number of 
Gaussians for each line are listed in Table [1] Note that we only use the spectral region of 
1450-1580A to fit the CiV line, as shown in Section EJl We don't measure Nv, Sill, Cll] 
and Mgll (also labelled in Figured]), because these lines are severely blended with either 
strong Lya or strong Fell bump. Due to the uncertainty in the origin of the I6OOA bump 
(see N06), it is hard to determine the line profile of Hell. We don't investigate Hell and 
Olll] too. 

The errors of the line fluxes shown in Figure E] are propagated from the uncertainties 
of the composites. Each composite is built from more than 2000 quasars, the flux errors are 
therefore very small. The errors include none from the uncertainties in continuum estimate, 
which is difficult to quantify. However, the continuum placement is important for weak 
line estimates. In the spectral region of 1690-2030A, several weak lines are mildly blended. 
To decompose these lines, we select spectral regions near 1695A and 2030A, where the 
emission line contribution to the total flux appears to be small, as continuum windows. 
We mask the region of 1760-1800A, where Fell emission is prominent, then simultaneously 
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fit the NIV, Nlll], SiII+[NeIIl] and ClII]+SiIII] lines. Figure [U shows an example for the 
spectral decomposition. Our best-fit model reproduces the composite. Ol+Sill in small BAI 
composites appears to lie on the red wing of the strong line complex Lya+Nv, we therefore 
adopt local continuums around the line as continuum windows (see Tabled]). Nevertheless, 
the high continuum around Ol+Sill at small BAI might be partly contributed by the broad 
profile of Ol+Sill. If this is the case, the Ol+Sill intensity at small BAI is underestimated 
and more detailed decomposition is required to study its dependence on BAI. 



3. BAI and Metallicity Correlation 

The line fluxes as a function of BAI are shown in the left panel of Figure [31 One can 
see the line intensities are correlated with BAI, and there seems to exist opposite trends 
between strong and weak lines. Strong lines (CiV, OVI+Ly/3, ClII]+SiIII] and SilV+OlV]) 
either decrease or remain constant with the increase of BAI, while weak lines (the other 
lines shown in Figure [3]) are dramatically enhanced at large BAI. The opposite trends also 
appears to exist in the lines, which are not measured: the strong line Mgll and the weak 
lines Sill, Si III] and Cll] (see Figured] and Wangll). It is well known that line equivalent 
widths decrease with quasar luminosity, known as the Baldwin effect (Baldwin 1977; Dietrich 
et al. 2002). If BAI was correlated with luminosity, one would expect that line intensities 
are correlated with BAI. We compute the median luminosity of each subsample, and find 
that quasars in the largest-BAI subsample are, on average, only 45% more luminous than 
those in the smallest-BAI subsample. Apparently, BAI is weakly related to luminosit}0 and 
Baldwin effect is not important for the trends reported in this paper. 

We then show the line ratios of other ions to C IV as a function of BAI in the right 
panel of Figure [3] All line ratios, except O VI+Ly/3 / C IV, positively correlate with BAI with 
a high statistical significance. Among these ratios, SiIV+OlV]/ClV, AllIl/CiV, Cll/CiV 
and ClII]+SiIII]/ClV have been demonstrated, by the photoionization model, to signifi- 
cantly increase with increasing metallicity (e.g. N06). This means that there exists strong 
correlation between BAI and metallicity. The intensities of the strong metal lines depend 
weakly on metallicity because of the thermostat effect introduced by the strong cooling lines 



^N06 found a significant correlation between CIV blueshift and luminosity. We note that they divided 
the sample into subsamples in a very different way from us. Most of quasars in their sample are assigned 
into one or two subsamples, while our subsamples are equal size. It might be the most important reason for 
this difference. 
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(Korista et al. 1998; Ferland 1999). In metal-rich gas, Civ and Ovn 

the most important 

coolants, even tend to decrease with increasing metallicity because their relative importance 
as coolants decreases. While, the weak lines can prominently increase with metallicity be- 
cause the number density of corresponding ions increases and the thermostat effect for them 
is trivial. Therefore, the metallicity-BAI correlation can well explain the trends of the line 
intensities shown in the left panel of Figure [31 We note that metallicity can not be the only 
factor affecting line emission, other factors, such as the ionizing continuum and the overall 
properties of BELR gas, also play an important role. 

Considerable efforts had been devoted to developing line ratios, such as Nv/CiV and 
SiIV-|-OlV]/ClV, as metal abundance diagnostics (Hamann et al. 2002; N06). Here, we 
adopt SiIV+OlV]/ClV as a metallicity indicator. Two theoretical line ratio-metallicity re- 
lationships, calculated by photoionization models using two types of ionizing continuums, 
are adopted (N06, Table 10). As shown in Figure HI the correlation between the inferred 
metallicity and BAI is weakly dependent on the ionizing continuum. When BAI increases 
from 0.49 to 0.76, the metallicity remarkably increases from 1.7 to 6.9^0. 

We use this indicator rather than that involving Nv because SilV+OlV] is a well iso- 
lated feature and the origin of the Nv line might be complex (Wang et al. 2010a). The 
theoretical calculation in N06 shows that Cll/CiV, AllIl/CiV and ClII]+SiIII]/ClV are 
also correlated with metallicity. For comparison, we also show the metallicities inferred from 
ClII]-|-SiIII]/ClV and AlIIl/CiV as a function of BAI in Figure H The measured ratios 
of ClI to Civ for the two smallest-BAI subsamples are less than the theoretical ratios at 
0.2Zq, which is the lowest metallicity adopted in the theoretical calculation, so we don't 
calculate the 'ClI metalliciy'. One can see that the 'Aim metallicity' is in good agree- 
ment with the 'SilV+OlV] metalliciy', implying that the AlIII measurement is reliable. The 
'ClII]+SiIII] metallicity' also increases with BAI, although it is apparently lower than the 
'SiIV-|-OlV] metalliciy'. N06 also found that 'ClI and ClII]+SiIII] metallicities' are much 
lower than 'AllII and SilV+OlV] metallicities'. ClI and ClII]+SiIII] are emitted by gas 
with different ionization degree or gas density from CiV. This indicates that Cll/CiV and 
ClII]-|-SiIII]/ClV are sensitive to the BELR properties adopted by photoionization model, 
such as the density and spatial distributions of BELR clouds (see NOG for more details). In 
contrast, SiIV-|-Olv]/ClV is less sensitive to the variation in the BELR structure, and thus 
is a more robust metallicity indicator. 



^Notc that the measured intensity of OVI+Ly/3 is severely affected by Lya forest, its intrinsic intensity 
should be much higher than the measured value. 
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4. Summary and Discussion 

We measure various line fluxes and tlie bluesliift and asymmetry index (BAI) of tlie 
Civ line from composite quasar spectra built from the Sloan Digital Sky Survey. Most 
of investigated line-flux ratios increase signiflcantly with BAI. Based on phototionization 
models, we interpret these dependencies in terms of stronger outflow in higher gas metallicity 
environment. This correlation can also account for the positive trends for weak lines and 
negative trends for strong lines, which are observed in the composites. We then use the 
flux ratio of SilV+OlV] to CiV as abundance diagnostics and flnd that metallicity increases 
more than 3 times as BAI increases from 0.49 to 0.76. 

We recall that a remarkably similar correlation has been well established in radiation- 
driven stellar wind (Massey 2003). High metallicity leads to large gas opacity, thus large 
absorption and acceleration. Recent studies on quasar outflows also favor the paradigm 
of radiative acceleration (Wangll). So the metallicity-outflow correlation can be naturally 
understood in terms of the metallicity-acceleration relationship. Another possibility is that 
starburst in the galactic center can produce abundant metals and drive the gases to fuel the 
black hole (Wang et al. 2010b). Because outflows are boosted at high accretion rates, the 
BAI— metallicity relationship can also be produced. The two mechanisms may both work. 
Although it is difficult to disentangle the predominant physical mechanism from the current 
data, the first mechanism is apparently more direct. Whichever mechanism dominates, star 
formation occurred in the host galaxies, which produces metals, has an important impact 
on the driving of quasar outfiow. 

The quasar metallicity may dramatically affect the host galaxy properties via feedback. 
On one hand, outflows from the quasar could exert enough force on the ISM to remove 
some of the material from the galaxy, consequently inhibiting the growth of the galaxy. This 
process would be more effective at higher metallicity according to our results. On the other 
hand, metal-rich outflows and ejected ISM may chemically enrich the IGM as the BELR 
metallicity is much higher than the solar value (Hamann et al. 2002). The enrichment 
makes the radiative cooling of the IGM more efficient, and further enhances star formation 
in galaxies and the growth of their central black holes, in contrast to the momentum feedback 
process that blows away the ISM. The enrichment rate is proportional to both the mass loss 
rate of outflows and metallicity, and therefore increases with metallicity in a manner stronger 
than a linear relationship. This feedback process is particularly important at high redshifts 
because 1) the intense star formation at high redshifts may form a substantially metal-rich 
environment in the galactic nuclei; 2) at high redshifts, IGM is very metal poor and its 
cooling is less effective. 
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Fig. 1. — The composite quasar spectra, in the order of increasing BAI, are plotted as red, 
orange, green, cyan and blue, respectively. Most of the emission features are labelled. The 
dashed lines indicate the rest frame wavelength of the corresponding emission lines. The 
inset shows a demonstration of spectral fitting in 1700A ~ 2300A region for the middle-BAI 
composite (green line). The black line shows the best-fit model, the purple dashed lines show 
the fitting results for each emission line. The model is described in detail in Section 12.21 
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Table 1: Continuum windows and number of Gaussians 



Lines 


A,o(A) 


A^.(A) 


Gaussians 


OVI+Ly/3 


1012 


1055 


1 


Ol+Sill 


1286 


1321 


2 


Cll 


1321 


1353 


2 


Siiv+Oiv] 


1353 


1450 


3 


CIV 


1450 


1695 


3 


NIV 


1695 


2030 


1 


Nlll] 


1690 


2030 


1 


Si 11+ [Ne III] 


1690 


2030 


1 


Aim 


1690 


2030 


1 


Ciii]+Siiii] 


1690 


2030 


2 



Note. — A/o and \hi indicate the center wavelengths of two continuum windows. 




Fig. 2. — Number of quasars in each wavelength bin of the composites as a function of 
wavelength. The color code is the same as those in Figure [TJ 
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Fig. 3. — Left panel: Various line flux as a function of BAI. Right panels: The flux ratio 
of other lines to C IV as a function of BAI. All these quantities are measured from the 
composites. Note that O VI+Ly/3 are severely affected by Lya forest. 
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Fig. 4. — Metallicity inferred from different indicators (as indicated in the panel) as a 
function of BAI. Please see Nagao et al. (2006) for the details about the theoretical line 
ratio-metallicity relationships for different UV continuum models. 



